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We present a mechanochemical model for myosin V, a two-headed processive motor protein. We 
derive the properties of a dimer from those of an individual head, which we model both with 
a 4-state cycle (detached, attached with ADP.Pi, attached with ADP and attached without nu- 
cleotide) and alternatively with a 5-state cycle (where the power stroke is not tightly coupled 
to the phosphate release). In each state the lever arm leaves the head at a different, but fixed, 
angle. The lever arm itself is described as an elastic rod. The chemical cycles of both heads are 
coordinated exclusively by the mechanical connection between the two lever arms. The model 
explains head coordination by showing that the lead head only binds to actin after the power 
stroke in the trail head and that it only undergoes its power stroke after the trail head unbinds 
from actin. Both models (4- and 5-state) reproduce the observed hand-over-hand motion and 
fit the measured force-velocity relations. The main difference between the two models concerns 
the load dependence of the run length, which is much weaker in the 5-state model. We show 
how systematic processivity measurement under varying conditions could be used to distinguish 
between both models and to determine the kinetic parameters. 



Introduction 

Myosin V is a motor prote in involved in differen t forms 
of intrace llular transport (|Reck-Peterson et all 1200ft 
IValeLl200^ . Because it was the first discovered processive 
motor from the myosin superfamily and due to its unique 
features, including a very long step size, it has drawn a lot 
of attention in recent years and now belongs to the best 
studied motor protein s. The experiments have character 
ized i t mechanically ijM ehta et all Il999t iPurcell et al 



iRief et all 120001: iRock et all 12003 IVeigel et al 
_ ._, . biochemically foe L a Cruz et all l2000ai 11999 . 
2000U IPurcell et alll2002tlYengo et all 12002ft. optically 



zuuub: irurcen et auiauuaiYengo et auizuuai. optically 
(lAli et all 120021 iForkev et all 12003 lYildiz et all 1200.1 



and s tructurally four gess et a ll l2002t ICoureux et al 
12003 IWalker et all l2000t IWang et all 12003 . These 
studies have shown that myosin V w alks along actin fil- 
aments in a hand-over-hand fashion (lYildiz et all 12003 
with an average step size of about 35 nm, roughly corre 



sponding to the periodicity of actin f ilaments (lAli et al 
2002 } iMehta et all 11993 IRief et all Eiofj IVeigel et al 



2002p . a stall force of around 2 pN llRief et " all 12000ft 



and a ru n length of a few micro ns llBaker et al 
IRief et all Efflfj ISakamoto et all r2 



120041 

Under physi- 
ological conditions, ADP release wa s shown to be the 
time limiting step in t he duty cycle foe La Cruz et all 
1999: IRief etail 12000). Two stages of the power stroke 
have been resolved: one about 20nm, possibly connected 
with the release of phosphate, and another one of 5nm , 
proba bly occurring upon release of ADP l|Veigel et all 
2002). Despite all this progress, the definite answer to 
the questions how the mechanical and the chemical cycle 
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are coupled and how the heads communicate with each 
other to coordinate their activity has not yet been found. 

Theoretical models for processive molecular mo- 
tors can follow different goals. What most models 
have in common is that they identify a few long- 
living states in the mechanochemical cycle and as- 
sume stochastic (Markovian) transitions between them. 
The differences between models start in the way these 
states are chosen. An approach that has been ap- 



plied to myosi 


n V (Kolomeiskv and Fisher, 2003 


), ki- 


nesin llPeskin and Osterl 




Schicf and Howard, 


2001: 


iThomas et all 


12002ft. as well as to other biological 



mechanisms of force generation , including actin poly- 
merization llPeskin et alj. ^93) and RNA polymerase 
l|Wang et all Il998ft . models the motors as stochastic 
steppers. These models describe the whole motor as an 
object that can go through a certain number of confor- 
mations (typically a few) with different positions along 
the track. After the completion of one cycle (which is, 
in models for myosin V and kinesin, tightly coupled to 
the hydrolysis of one ATP molecule), the motor advances 
by one step. All steps are reversible and at loads above 
the stall the motor is supposed to walk backwards and 
thereby regenerate ATP. The approach has been partic- 
ularly useful for interpreting the measured force- velocity 
relations and relating th em to the kinetic parameters 



and positions of substeps IIFis 


ter and KolomeiskyM200ll: 


Kolomeiskv and Fisherll2003 


Schief and Howard! 12001ft. 



A limitation of such models is that they assume coordi- 
nated activity of both heads rather than explaining it. 
They also assume that the motor strictly follows the reg- 
ular cycle and there is no place for events like steps of 
variable length and dissociation from the track, although 
the latter can be incorporated into the models by propos- 
ing a different dissociation rate for each state in the cycle. 

In this Article we present a physical model for the pro- 
cessive motility of myosin V. The basic building block of 
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our model is an individual head, which we model in a sim- 
ilar w ay as the models for conventional myosins do llHill 
albeit with different rate constants. The head is 
connected to the lever arm, which we model as an elastic 
rod, wh ose geometry we infer from electron micr oscopy 
studies {Burgess et all 120021 IWalker et all l200(t . The 
two lever arms are connected through a flexible joint and 
this is the exclusive way of communication between them. 
We will derive the properties of the dimer from those of 
the individual head. 



The Model 

To describe each myosin V head we use a model based 
on the 4-state cycle as postulated bv iLvmn and Tavlorl 
ill 97 if) and used in many quantitative muscle models 
l|Hillt Il974j) (Fig.EJV). We restrict ourselves to the long- 
living states in the cycle: detached with ADP.Pi, bound 
with ADP.Pi, bound with ADP, detached with ADP and 
bound without a nucleotide. The bound state with ATP 
and the free stat e with ATP have both be en found to be 
very short-lived l)De La Cruz et all Qj)99) and we there- 
fore omit them in our description, i.e., we assume that 
binding of ATP to a bound head leads to immediate de- 
tachment and ATP hydrolysis. The detached state with- 
out a nucleotide is very unlikely to be occupied because 
of the low transition rates leading to it and we omit it 
from our scheme as well. 

One question that has not yet been definitely answered, 
is whether Pi release occurs before or during the power 
stroke, i.e., whether a head which is mechanically re- 
strained form conducting its power stroke can release Pi 
or not. The 4-state model assumes a tight linkage be- 
tween the Pi release and the power stroke. While the 
4-sta te model has b een successfully applied to myosin II 
(e.g., iDukd ||1999|) : IVilfan and Dukel l|2003jn . recent ex- 
perimental evidence suggests th at the lead head can re- 
lease Pi before the power stroke ( Ros enfeld and Sweenevl 
l2004j) . We therefore also discuss an alternative 5-state 
model. In the 5-state model we introduce an addi- 
tional state ADP' in which the phosphate is already re- 
leased, but the lever-arm is still in the pre-powerstroke 
state. The next transition, ADP release, however, is still 
linked to the completion of the full power-stroke. This 
is necessary in order to explain head coordination and 
also in agreement with experiments that show a strain- 
dependenc e in the ADP rele ase rate in single-headed 
molecules (jVeieel et all l2002j) . The extended duty cy- 
cle of a head is shown in Fig. 03. 

A head always binds to an actin subunit in the same 
relative position. In each state, the proximal end of the 
lever arm leaves the head in a fixed direction in space, 
determined by the polar angle </> towards the filament 
plus end and the azimuthal angle = 9$i of the actin 
subunit i to which the head is bound. The geometry of 
the molecule and the angles were inferred from images 
obtained with electron microscopy {Burgess et all l2002t 
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FIG. 1 The myosin V dimer is modeled as two heads, each 
connected to a lever arm which leaves the head at a certain 
angle <f>, depending on the state of the head. The two lever 
arms, modeled as elastic beams, are connected with a flexible 
joint, which is also connected to the external load. 



+ATP 




FIG. 2 A) The mechanochemical cycle of each individual 
head. The head attaches to actin in the state with ADP and 
Pi bound on it, undergoes a large conformational change upon 
Pi release, another smaller conformational change upon ADP 
release, then binds ATP and enters the very weakly bound 
state, which dissociates quickly. B) The mechanochemical 
cycle in the 5-state model. In this scenario, the phosphate 
release and the power stroke are two separate transitions. 
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TABLE I Geometric parameters of a myosin V head (see also 
Fig- Q for their definition). 



Lever arm length L 



26 nm 



Lever arm start R 8 nm 

Lever arm start 5adp.pi Onm 

Lever arm start ^ADP.apo 3.5 nm 

Angle ADP.Pi 0adp.pi 115° 

Angle ADP 0adp 50° 

Angle apo apo 40° 



I Walker et all 120(1(1) . They are summarized in Table [I] 
In our calculations we assume a 13/6 periodicity of the 
actin helix (6 rotations per 13 subunits), which means 
O = 2tt x 6/13. 

We assume that the lever arm has the properties of a 
linear, uniform and isotropic elastic rod, described with 
the bending modulus EI. Then the local curvature k is 
determined from M = EIk, where M is the local bend- 
ing moment (torque). The lever arms from both heads 
are joined together (and to the tail) with a flexible joint 
which allows free rotation in all directions. For a certain 
configuration of chemical states, binding sites of both 
heads and a given external force, the three-dimensional 
shape and the bending energy of both lever arms can 
be calculated numerically as described in the Appendix. 
Some of the calculated shapes are shown in Fig. 

We calculate the free energy of a dimer state as 



where AU denotes the change in elastic energy of the 
dimer and FAx the work performed against the external 
load. 

The second exact statement can be derived by mul- 
tiplying together the detailed balance conditions for a 
monomer in the absence of any external force along a 
closed pathway in Fig. After one cycle the free en- 
ergy of the bound monomeric head returns to its initial 
value, while the total free energy change in the system 
equals the amount gained from the hydrolysis of one ATP 
molecule. The resulting relation reads 

fc+ATP [ATP] 



k- a fc+pi [Pi] ^+adp [ADP] k°_ ATp 

^GATP AG° [ATP1 

= e k B T — e k B T — - — (3) 

[ADP] [Pi] 1 ' 

and provides an important constraint on the kinetic rates 
of the model. In the 5-state model, we obtain an equiv- 
alent equation, 



fc+ATP [ATP] 
k- a fc+p; [Pi] k°_ PS fe + ADP [ADP] k°_ ATp 



AG 
— & h r T . 



[ATP] 
[ADP] [Pi] 



(4) 



A similar statement also holds for the rates along the 
inner loop in the reaction scheme, which involves attach- 
ment, power stroke and detachment, all in the ADP state. 
Because we assume that the detachment rate in the pre- 
powerstroke and the post-powerstroke state are both the 



G = G 1 + G 2 + U 1 + U 2 + Fx 



(1) same (k'_ A ), the relation reads 



where G\ and G 2 are the intrinsic free energies of both 
heads (which depend on the chemical state of the head 
and the concentrations of nucleotides) , U% and U 2 are the 
energies stored in the elastic deformation of each lever 
arm, and Fx is the work done against the external load 
(x denotes the coordinate of the flexible joint along the 
filament axis with positive values towards the plus end, 
while positive values of F denote a force pulling towards 
the minus end, against the direction of motion of an un- 
loaded motor). 



Transition rates 

There are two exact statements we can make about 
the kinetic rates of the duty cycle that follow from the 
principle of detailed balance. The first statement relates 
the forward and the backward rate of any reaction to the 
free energy difference between the initial and the final 
state. For any transition the principle of detailed balance 
states that 



k+j 

k—i 



k° ■ ' 



AU+FAx 



(2) 



If " k° 

1.0/ 1.0 
re +A -PS 



(5) 



When it comes to the actual force dependence of tran- 
sition rates we have to rely on approximations. An ap- 
proach that is most widely used when modeling motor 
proteins, but also other conformational changes, like the 
gating of ion c hannels, inv olves the Arrhenius theory of 
reaction rates l)Hill Il974|) . It proposes that the pro- 
tein has to reach an activation point (x a ) somewhere be- 
tween the initial (xi) and the final state (xf) by thermal 
diffusion, but completes the reaction rapidly after that. 
Therefore, the force dependence of the forward rate can 
be modeled as 



e- 



k°_ ie - 



(6) 



where U (x) means the total potential (bending of both 
lever-arms and work done against the external load) 
which a head has to overcome to bring the lever arm 
angle into a given state. We use the variable e to denote 
the relative position of the activation point between the 
initial and the final state, so that x a = (1 — e)x t + exf. 
Unless otherwise noted, we will assume e = 0.5. Not 
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FIG. 3 Calculated shapes and bending energies of dimers, bound i subunits apart [i — —2,2, . . . , 15) and in different states: 
first in post-, second in the pre-powerstroke state (upper row), both in the post-powerstroke state (middle row) and both in the 
pre-powerstroke state (bottom row). Each configuration is shown in side and front view. If both heads are in the same state 
(bottom two rows) there is a significant cost in elastic energy needed to buckle one of the lever arms. Binding of the lead head 
before the trail head undergoes the power stroke is therefore unlikely. 



precisely identical, but useful for practical purposes is 
also the approximation U(x a ) = (1 — e)U(xi) + eU(xf). 
Therefore we get the following expression for the force- 
dependence of the transition rate: 



k 



+i 



k° 



(7) 



For reactions that involve the binding and unbinding 
of a head, Eq. [2] is valid, but one expects the activa- 
tion point to be much closer to the bound state. The 
strain-dependence of the detachment rate for heads in the 
ADP and ATP.Pi state has not yet been measured and 
we therefore neglect it, assuming that the detachment 



rate is force-independent, fc_ 



-A- 



The attachment 



rate then relates to the potential difference as 



'.,0 -7 



(8) 



Choice of kinetic parameters 

Some of the transition rates in the cycle are well 
known from the literature. &_adp is the limiting 
rate both for running myosin V molecules and for 
single-headed constructs at low ATP concentrat ions. 
The measured valu es are 13 s- 1 llRief 'et all 12(10(1) for 
dimers and 12 s" 1 jDe La Cruz et all Il999t), 13-22 s" 1 
llTrvbus et ail 11999(1 . and 4.5-7 s 1 <|Moilov and Veigel 
2003) for monomers. Because the actual rate in a 
dimer is slowed down as compared to the monomer, 
we use the value k°_ ADP = 20 s -1 . The reverse rate, 
fc+ADP can be determined from the inhibitory effect 
of ADP on th e velocity and has been estimated as 
n.GuM^s- 1 llDe La Cruz et all Il999l). ^SuM 'V 1 
llRief et all 1200(1) . MuM^s" 1 iWang et all 1200(1) . 



Equally well 
binding, fc+ATF 



0.9 /x rvrV 1 



known is the rate for ATP 
which has b een me asured as 
|5e La Cruz et all Il999t iRief et all 
1200(1) . 0.6-1 .5 //.TVr's- 1 l|Veigel et all 12002^ For 
the Pi release rate the estim a tes range from 
fc-Pi > 250s" 1 tDe La Cruz et all Il999|) to 110s" 1 
l|Yengo and Sweenevl . 12004^ . We therefore use the value 
fc-Pi = 200 s" 1 . 

There is some more discrepancy between the cur- 
rent values for the release rate from actin in 
the ADP state. While direct measu rements gave 
fc'_ A = 0.032 s- 1 llDe La Cruz et all Il999f) and 0.08 s" 1 
l|Yengo and Sweenevl 120041) . a recent estim ate from the 
run l ength led to a higher value of 1.1s -1 ((Baker et all 
I2004D . We use an intermediate value of fcl A = 0.1s . 
For the attachment rate in the ADP s tate, we set fc+ A sa 
fc^ A , based on kinetic measurements l)De La Cruz et all 
1999). 

This leaves us with a total of 4 unknown kinetic rates, 
of which 3 need to be estimated from the measured step- 
ping behavior and run length data, while one can be de- 
termined from Eq. |3] 



Results 

Choice of the value for the bending modulus 

There are two ways to estimate the bending stiffness of 
the myosin V lever arm - one from its structure and anal- 
ogy with similar molecules and the other one from the ob- 
served behavior of the dimeric molecule. The lever arm 
consists of 6 IQ motifs, forming an q-helix, surrounded 
by 6 calmodulin or other light chains ijTerrak et all 120031 
IWang et all 12003). One possible estimate for the stiff- 
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TABLE II Kinetic parameters of the model 



Parameter Value Source 

4-state 5-state 



fc+A actin binding with ADP.Pi 
k-A actin release with ADP.Pi 
actin binding with ADP 
actin release with ADP 

Pi release 



K +A 



e-Pi 



ft +PS 



activation point 
Pi binding 
power stroke 
reverse stroke 



k°_ ADP ADP release 

fc +ADP ADP binding 

fc+ATP ATP binding, actin release 

k°_ ATP actin binding with ATP release 



5000 s" 1 


5000 s" 


Is" 1 


50 s" 1 


5000 s- 1 


5000 s" 


0.1 s" 1 


0.1s" 1 


200s -1 


200 s" 1 



est. from run length 
est. from run length 



0.3 

10" 4 /iM" 



20 s~ 



12 fiMr 



0.7/dVT 



0.07 s- 



-i 



10" 
10 4 s" 
0.05 s 
20 s" 1 



12 /dVT 



0.7/iM" 



k' 



+A 



ijPe La Cruz et all 



0.032 s" 1 JPe La Cruz et"ail 
1.1s" 1 llBaker et all l2004fl 



l|De_I^g L j2nrz_ej i _all 

110 s" 1 iYengo and Swccnevl 



> 



250 s" 



1199' 

228 s" 1 jRosenfeld and Sweenevt kOO^ ^ 
F-v relation at high loads 
guess 
guess 

fc°pg/fc° ps from the stall force 
fc-AD P = 13 s" 1 for dimers jRief et all 
2000) 

12.6 /xM~ 1 s" 1 
14 /iM-V 1 



1.2 s" 



JPe La Cruz et all U '.)!): 
iWang et alll2000l) 

0.9 ^M" 1 ^ " 1 JPe La Cruz et all |199 
iRief et all |2000D 0.6 - 1.5 ^M" 1 
llVeieel et alll2002l) 
Eq.EO Eq.H 



ness of the lever arm can be obtained by approximat- 
ing; it with a coiled- c oil do main, as has been done by 
iHoward and Spudicbl 1^)96?). Generally, the stiffness of a 
seminexible molecule is related to its persistence length £ p 
as EI = £ p fcsT '. Howard and Spudich estimated the per- 
sistence length of a coiled-coil domain as 100 nm, which 
yields EI « 400pNnm 2 . Other researchers report val- 
ues of £„ — 130 nm for myos in llHvidt et all 11982 ) and 



150 nm for tropomyosin ( Phillip s and Chackol 



1996; 



ISwenson and StellwagenL 11989^) . 

On the other hand, we can estimate the stiffness from 
the force a lever arm has to bear under conditions close 
to stall. We do this by calculating the distribution of 
binding probabilities to different sites at F = 1.8 pN, 
which is close to stall force. We assume that the binding 
rate to each site is proportional to its Boltzmann weight, 
exp(— G/fcflT), which is equivalent to assuming that the 
activation point of the binding process is close to the final 
state and that the reverse reaction (detachment in the 
state with ADP.Pi) has no force-dependence in its rate. 
The expectation value of the binding position of the lead 
head relative to the trail head is shown in Fig.0] It shows 
that a stiffness of EI > 1000 pN nm 2 is necessary to allow 
stepping at loads of this magnitude. 

For these reasons, we use the value EI = 1500 pN nm 2 . 
This corresponds to an elastic constant (measured at the 
joint) of 




2000 



EI (pN nm ) 



k = 3EI/L 3 = 0.25pN/nm 



(9) 



FIG. 4 The average step size under a load of F = 1.8 pN 
as a function of the lever arm elasticity EI. The step size 
was calculated from attachment probabilities of the lead head 
(APP.Pi state) relative to the bound trail head (APP state). 



The elastic constant for longitudinal forces (with respect 
to the lever arm) is much higher. If we approximate 
the lever arm with a homogeneous cylinder of radius 
r = 1 nm, we can estimate it as = 4EI/(r 2 L) — 
230pN/nm. We therefore neglect the longitudinal exten- 
sibility of the lever arm in all calculations. 

A similar value (EI = 1300 pNnm 2 ) has also been ob- 
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tained by analyzing data from optical trap experiments 
on single-he aded myosin V mole cules with different lever 
arm lengths l|Moore et al.Ll206"Ih . Even though it is some- 
what large r (about 3 times) than the va lues estimated for 
myosin II ijHoward and Spudichl . Il99r^ . there is no solid 
evidence that the structures with different light chains 
have the same bending stiffness. On the other hand, 
there could have been some evolutionary pressure to in- 
crease the lever arm stiffness, as it is directly related to 
the stall force of myosin V. While we are not able to give 
a definite answer to the question whether the lever arm 
behaves like a uniform elastic rod or whether there is a 
pliant region close to the head, we favor the first hypoth- 
esis because the estimated lever arm elasticity already is 
more than sufficient to explain the mechanical properties 
of the dimeric molecule. 



Step size distribution 

Figure [3] shows the energies stored in the elastic dis- 
tortions of the lever arms of both heads in the pre- 
powerstroke or the post-powerstroke state. For example, 
if the first head is in the ADP.Pi state and the second 
head binds before the first one undergoes a power stroke, 
this is connected with an energy cost of 6.6 fcsT. The 
attachment rate of the lead head before the power stroke 
in the trail head is therefore more than 100 times slower 
than after the power stroke. 

Because the lead head normally attaches to actin while 
the trail head is in the ADP state, we can determine the 
probability that the lead head binds to an actin site i 
subunits in front of the trail head from the Boltzmann 
factors formed from the bending energy in the final con- 
figuration, Pi oc exp(-([/i + U 2 )/k B T). Here £7i + U 2 
denotes the sum of elastic energies stored in both lever 
arms if the trail head is in the ADP state and the lead 
head in the ADP.Pi state, bound i sites in front of the 
trail head. The resulting distributions for different lever 
arm lengths are shown in Fig. For the lever arm con- 
sisting of 6 IQ motifs, the result is a mixture of 11 and 13 
subunit steps, whereby 13 subunits dominate. Azimuthal 
distortion plays a major role in the bending energy, there- 
fore binding is only likely to sites 2, 11, 13 and 15, on 
which the azimuthal angles of both heads differ by not 
more than 27°. 



The gated step in the cycle 

A question that has been a subject of intense discussion 
is which step in the cycle is deciding for the coordination 
of the two heads. A currently often favored hypothesis 
proposes that the lead head undergoes its power stroke 
immediately after binding, thereby storing energy into 
elastic deformation of its lever arm and releasing it after 
the unbinding of the trail head. An alternative hypothe- 
sis proposes that the release of the rear head is necessary 
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Step size (actin subunits) 

FIG. 5 Step size distribution for 4 different lever arm lengths 
L: lOnm (2IQ), 18nm (4IQ), 26nm (6IQ) and 34nm (8IQ) and 
no external load. The histograms show the probability that a 
lead head (ADP.Pi state) will bind i sites in front of the trail 
head in the post-powerstroke ADP state. The probabilities 
were determined from the Boltzmann factors, resulting from 
the elastic distortion energy of the configuration. Azimuthal 
distortion plays a crucial role role in determining the step size, 
which is the reason why the binding is always concentrated on 
sites 2, 11, 13 and 15. Taking into account the fluctuations 
in the actin would lead to a broader distr i butio n, in better 
agreement with experiments llWalker et all 12000). 

for the power stroke in the front head. As we will show 
below, our model favors this picture. In the 4-state sce- 
nario, this implies that the lead head is waiting in the 
ADP.Pi. In the 5-state scenario it is in the ADP' state 
(the pre-powerstroke ADP state). The trail head spends 
most of its cycle in the ADP state in both scenarios at 
saturating ADP concentrations. 

Because this model challenges the currently prevail- 
ing view, we should first critically review the arguments 
supporting it. One argument includes the direct observa- 
tion of telemark-shaped molecules, with the leading head 
leaning for ward and then the lever arm tilted strongly 
backwards ( Walker et~al l l200fl) . A more detailed image 
analysis, however, showed that the conver ter of the lead- 
ing h ead is in the pre-powerstroke state ijBurgess et all 
2002) . Another p iece o f evidence comes from experiments 
bv iForkev et alJ l)2003|) which show a fraction of tags on 
the lever arm (30-50%) that do not tilt while moving, 
but again the data provide no conclusive proof because 
the method does not allow detection of tilts symmetric 
with respect to the vertical axis. To conclude, one cannot 
say that the present experimental evidence excludes any 
of the two hypotheses about the moment of phosphate 
release and of the power stroke. 

From the theoretical side, we will argue that in a 
model with linear elasticity the mechanism with imme- 
diate power stroke in the lead head cannot work under 
loads for which the motor is known to be operational. 
It is known that the monomeri c constructs of myosin V 
undergo a normal duty cycle l)De La Cruz et all Il999t 
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lYengo et jH EOOl, 

which means that no step in the cycle 
requires mechanical work from the outside for its comple- 
tion (which would be, for example, the case if the head 
needed to be pulled away from actin to complete the cy- 
cle). This excludes the possibility that the free energy 
gain connected with binding and the power stroke ex- 
ceeds AGatp = lOOpNnm, the total available energy 
for one cycle. Because this and other transitions in the 
cycle need to be forward-running, we use the still con- 
servative estimate that the free energy gain from binding 
and the power stroke cannot exceed 80pNnm. On the 
other hand, we can estimate the free energy that would 
be necessary for a head to bind to a site 13 units ahead 
and then undergo a conformational change. The amount 
of energy needed to bring the dimer into the hypothetical 
state with both heads in the post-powerstroke state and 
a strong distortion, especially of the leading lever arm, is 
plotted in Fig. The calculation shows that the bind- 
ing of the front head with the subsequent power stroke 
before the rear head detaches (for a load of F = 1.8 pN) 
is only possible for values of EI < 450pNnm 2 , which 
is inconsistent with the lower estimate based on the ob- 
served step size (Fig. 0J. Of course, we cannot rule out 
that there is some additional state in the middle of the 
power stroke which is occupied immediately while the 
lead head waits for the trail head to detach. But within 
the scope of the geometrical model with a single power 
stroke connected with the Pi release, we consider the sce- 
nario where the lead head instantaneously undergoes the 
power stroke without waiting for the detachment of the 
trail head unrealistic. 
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FIG. 6 The amount of energy needed for the binding of the 
lead head and the subsequent power-stroke, plotted against 
the lever arm elasticity. The load pulling on the tail is 
F = 1.8 pN. The lower curve shows the energy needed to 
pull the external load and distort the lever arms in order to 
bind the new lead head 13 sites in front of the trailing head. 
Note that most of this work will be performed by Brownian 
motion, but the potential well in the bound state still has 
to be strong enough to stabilize the bound state. The middle 
curve shows the energy needed mainly for the distortion of the 
lever arms when the lead head undergoes a power-stroke be- 
fore the trailing head detaches. Since the sum of both cannot 
be higher than 80 pN nm, we estimate that this hypothetical 
scenario would only be possible if the lever arm stiffness was 
EI < 450pNnm 2 . This is inconsistent with other require- 
ments of the model, so we rule this scenario out. 



Hidden power strokes in the dimer configuration 

An immediate consequence of the elastic lever arm 
model is that the tail position is mainly determined by 
the geometry of the triangle and less by the conforma- 
tions of individual heads. For a monomeric head or a 
dimer bound by a single head, the power-stroke upon 
ADP release has an x-component (in the direction of the 
actin filament) of about 3.3 nm (Fig. 01. If the lead head 
is attached, however, the power stroke as measured on 
the tail is reduced by about a factor of 50. The tail 
movement is also closely related to the force-dependence 
of transition rates, which means that transitions between 
states with both heads bound do not show any significant 
load dependence. In the kinetic scheme we use here this 
implies that the rates of ADP release and ATP binding 
(the two rate limiting steps at low or forward loads) are 
both con stant, in agreement with the flat F-v curve mea- 
sured bv lMehta etall l)l999|) . 



Force-velocity and run length curves 

The bending energies, calculated for each possible 
dimer configuration, and the transition rates were fed 




FIG. 7 For a single head, the ^-component of the power- 
stroke upon ADP release equals 3.3nm (for zero load). In 
the dimer with both heads bound, only 0.07nm of that power 
stroke reach the load. As a consequence, the load-dependence 
of transition rates between states with both heads bound is 
negligible. 



into a kinetic simulation to determine the average ve- 
locity of a dimeric motor and its dissociation rate from 
actin. The most probable kinetic pathway of the dimer 
is indicated by thick arrows in Fig. |SJ while the thin ar- 
rows indicate some of the possible side branches that can 
lead to dissociation. Figure ED shows the resulting force- 
velocity curves and Fig. ^2 the dissociation rates. 

An analytical solution of the 4-state model would, in 
theory, require solving the occupation probabilities for a 
system with about 6 + 8x3x3 = 78 states (6 states with 
one head bound, plus configurations with both heads 
bound, where each head can occupy 3 different states and 
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FIG. 8 Most probable kinetic pathways for a dimer in the 4-state (A) and in the 5-state model (B). The thick arrows denote 
the regular pathway and the thin arrows side branches that can result in dissociation from actin. Note that the simulation was 
not restricted to the pathways shown here, but included all possible combinations of transitions between monomer states. 



the relative positions of both heads can have 8 different 
values). Such a system could easily be solved numeri- 
cally, but would be too complex for obtaining an insight- 
ful analytical expression. However, we will show that a 
simplified pathway can already lead to expressions that 
agree reasonably well with simulation data and are there- 
fore useful for fitting model parameters to experimental 
data. 

In the following, we give approximate expressions for 
the most significant steps in the mechanochemical cycle. 
The average time it takes for a head in the state to 
bind an ATP molecule can be estimated as 

((+aTp) = ^^ 1 + ^±^2m (10) 

fc+ATp[ATPJ V fc_ADP / 

where the second term takes into account a reduction of 
the forward rate due to ADP rebinding. The second rate 
limiting process (especially at hight loads) is the release 
of phosphate. The average dwell time in the state with 
one head free and the other one in the ADP.Pi state is 



(*-Pl> 



1 



fc_ 



Pi 



(11) 



The third rate limiting step is the ADP release, with the 
time constant 



(t-ADp) 



1 



-ADP 



(12) 



With these three average dwell times, the motor velocity 
can be calculated as 



(d) 



(t-Pi) + (i-ADp) + (£+ATp) 



(13) 



where (d) denotes the average step size, which 
is about 35 nm. The individual rates that ap- 
pear in this expression can be estimated as 
follows: k—pi ~ fc° Pi exp(— Fe^pidps/fcsT) 
with dps = L(cos0adp — cos </>adp.pO + 5 and 

ADP « k°_ ADP exp(-AC/ _ADp/2fcsT) w 0.65/c°_ ADP . 
The results for two different ATP concentrations are 
shown in Fig. |5J^ and compared with a simulation 
result. The analytical expression reproduces the simu- 
lation result well, with a small deviation being mainly 
the result of alternative pathways, neglected force- 
dependence of the ADP release rate and variation in the 
step siz e. The expe ri menta l ly measured fo r ce-vel ocity 
curves l|Mehta et all Il999t lUemura et all 12004(1 are 
also well reproduced, although the experiments show 
a more abrupt drop in velocity at high loads, with no 
measurable effect up to about lpN. 

In the 5-state model the power-stroke can be fast 
and reversible, in which case the pre- and the post- 
powerstroke state can reach an equilibrium and the lim- 
iting rate is proportional to the probability of the post- 
powerstroke state 1/(1 + exp(F<ips/fc_BT)) - a signifi- 
cantly sharper load dependence than the 4-state model 



9 




FIG. 9 A) Force-velocity curves in the 4-state model, ob- 
tained from a stochastic simulation. The solid curve shows 
the values for lOOO^M ATP and the dashed curve for IfM 
ATP. Both curves are compared with the prediction of the 
simplified analytical expression (Eq. 1131 , dotted lines. The 
minor deviation is mainly due to cycles taking other pathways, 
neglected force-dependence of the ADP release rate and vari- 
ation in the step size. Note that the velocities above ~ 2.5 pN 
are not well defined because the dissociation time becomes 
comparable with the step time. B) Inhibition by ADP and 
Pi. The force-velocity relation with ImM ATP is shown by 
the continuous line. The dashed line shows the same relation 
with additional 10/iM ADP and the dotted line with ImM 
phosphate. The velocity reduction through ADP occurs at 
low or negative loads, while the inhibition by Pi only becomes 
significant close to stall conditions. 

(Fig.EJ. 



Inhibition by ADP and phosphate 

It is a well established observation that ADP can slow 
down myosin V by binding to heads in the state with no 
nucleotide and thereby preventing them from accepting 
an ATP molecule. The rate of ADP rebinding is already 
taken into account in the kinetic constants and the model 
naturally reproduces the observed behavior, as shown in 




F(pN) 



FIG. 10 Force-velocity relation of the 5-state model with 
ImM ATP (solid), ImM ATP+ lO/AI ADP (dashed) and 
1/iM ATP (dotted). Note the sharper drop at high loads as 
compared to the 4-state coupled model (Fig. |UJ . 

Figs. and ^] for the 4-state model and in Fig. f° r 
the 5-state model. Not yet experimentally investigated 
has been the inhibition by phosphate. Its intensity de- 
pends on the reverse power-stroke rate, which is one of 
the open parameters of our model. In the 4-state model, 
Pi re-binding is necessary for the reverse power stroke 
and therefore some inhibition effect can be expected at 
high loads. The simulation shows clearly that the phos- 
phate concentration has no effect on zero-load velocity, 
but it does slow down the motor close to stall (Fig. 03). 
A similar effect of Pi on isometric force h as also been ob- 
served in muscle llCooke and Patd. Il985h . In the 5-state 
model Pi rebinding is not mechanically sensitive and its 
effect is roughly force-independent. However, with the 
parameters chosen here, it is negligible. 



Three dissociation pathways 

As we can see from the kinetic scheme (Fig. |SJ, there 
are three significant pathways in the cycle that can lead 
to the dissociation of the myosin V dimer from an actin 
filament. The first pathway leaves the cycle if a dimer 
bound with one head in the ADP.Pi state detaches be- 
fore the second head can attach. The second pathway 
runs through a state in which the bound head releases 
ADP and binds a new ATP molecule before the free 
head can bind. With the third pathway we denote all 
processes that involve the detachment of a head in the 
A DP state. Th i s is t he pathway favored by recent results 
of iBaker et all ll2004lh Figure 1111 shows the dissociation 
rate, separated by contributions of the three pathways. 
They have the following characteristics: 

Pathway 1: With this pathway we denote the dissocia- 
tion of a head in the ADP.Pi state. Because this state is 
long-lived at high loads in the 4-state, but short-lived in 
the 5-state model, the resulting force-dependence of the 
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dissociation rate differs significantly in both scenarios. 
In the 4-state model, the contribution to the dissociation 
probability per step shows a strong load-dependence, but 
no significant dependence on the ATP concentration. It 
can be estimated as 



fc-A 
t_ Pi 



(14) 



with dps = L(cos0adp — cos </>adp.p0 + 5. The dissocia- 
tion rate is higher for positive loads. From the estimated 
run length at 1 pN load and saturating ATP concentra- 
tion of about 15 steps l)Clemen et all l2003|) . we can es- 
timate the unbinding rate as fc_A ~ Is -1 . In order to 
account for reported run lengths of over 50 steps at low 
loads, we tentatively assign fc° A w 5000 s -1 . 

In the 5-state model, the situation is reversed. There 
the dissociation process on path 1 takes place if the trail 
head releases ADP before the lead head releases Pi, which 
can happen in two different ways: on one the rate is 
approximately force- independent, on the other it grows 
with negative (forward) loads. In order to obtain a signif- 
icant contribution to the detachment rate on this path- 
way, we choose a higher detachment rate k~A than in the 

4- state model (50s _1 instead of ls^ 1 ). 

Pathway 2: Because the process of unbinding requires 
an ATP molecule, the per-step dissociation rate grows 
with the ATP concentration. In addition, it is propor- 
tional to the ratio of the ADP dissociation rate and the 
actin binding rate, A-adpA+a, which is higher for neg- 
ative (forward) loads. This holds in both the 4- and the 

5- state scenario. 

Pathway 3: The dissociation probability on pathway 3 
is proportional to the detachment rate in the ADP state, 
k'_ A . Of all three pathways, this one shows the weakest 
load-dependence, although it is higher for forward loads. 

We expect that systematic data on mean run length as 
a function of load and nucleotide concentrations will be 
helpful to determine the remaining model parameters. 
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FIG. 11 Dissociation rate of myosin V dimers from actin un- 
der a high (top) and a low (bottom) ATP concentration (4- 
state model). The continuous line shows the total dissociation 
rate, the dashed line the dissociation via pathway 1, the dot- 
dashed line via pathway 2 and the dotted line via pathway 
3. 
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Reverse stepping in the 5-state model 

As a consequence of both the reversibility of the power 
stroke and the slower dissociation rate at high loads, the 
motor can step backwards under loads exceeding the stall 
force (Fig. I14fl . Note that these steps are not the sim- 
ple reversal of forward steps (which would involve ATP 
synthesis) , but rather indicate a different pathway in the 
kinetic scheme, in which both heads stay in the ADP 
state and alternately release actin at the leading position 
and rebind at the trailing. The time scale of reverse step- 
ping is determined by the dissociation rate of a head in 
the ADP state, k'_ A , which we chose as 0.1 s _1 . With 
a higher value of k'_ A , especially for the pre-powerstroke 
state (so far we assumed that the rate is equal in both 
ADP states), faster stepping would also be possible, al- 
though there is an upper limit on k'_ A , imposed by the 
dissociation rate on pathway 3. 
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FIG. 12 Velocity (continuous, left scale) and mean run length 
(dashed, right scale) as a function of ADP concentration in 
the 4-state model for zero load and ImM ATP. 
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FIG. 13 Force- dependence of the dissociation rate in the 5- 
state model. The load dependence for positive loads is much 
weaker than in the 4-state model (Fig. li lt 
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FIG. 14 Reverse stepping in the 5-state model under a high 
load (4.5pN), 10 fiM ATP and 1 fiM ADP. There is also some 
creeping motion between the steps, which results from the 
attachment and detachment of the two heads on neighboring 
sites, and only takes place if myosin V is allowed to follow 
a helical path on actin. If binding is constrained to one side 
of the actin filament (like on a coverslip), then only regular 
reverse steps with the periodicity of the helix are observed 
(not shown). 



Discussion 

We used the geometrical data of the myosin V molecule 
as obtained from EM images to calculate the confor- 
mations and elastic energies in all dimer configurations. 
These data were first used in a model with a four-state 
cycle and subsequently in a five-state model. 

The first result, which follows directly from the bend- 
ing potentials and is independent of the underlying cycle 
is that the elastic lever arm model explains two key com- 
ponents of the coordination between heads: why the lead 
head does not bind to actin before the power stroke in 



the trail head and why it does not undergo its power 
stroke before the trail head detaches. It also allows us 
to calculate the distribution of step sizes. The results for 
different lever-arm lengths (Fig. |3J give realistic values, 
in agreement with step size and hel icity measurements 
(|AU et all l2002t iPurcell et all l2002(h even though they 
have a slight tendency towards underestimation and also 
show a narrower distribut ion than direct electron mi- 
croscopy observations ijWalker et al l l2000() . A possible 
explanation for the broader distribution than predicted 
by the model lies in the fact that in reality the actin struc- 
ture does not follow the perfect helix, as assumed in our 
mode l, but has angular dev iations of up to 10° per sub- 
unit ( Egel man et all ^982). Taking these fluctuations 
into account would clearly broaden the distribution of 
our step sizes, but alone it cannot explain the tendency 
towards longer steps. The most straightforward expla- 
nation for the longer steps is that the power stroke has 
an additional right-handed azimuthal component. Then 
the configuration with the lowest energy is reached if the 
lead head is twisted to the right relatively to the trail 
head, which is the case if it is bound further away along 
the helix. The observation that the actin repeat is often 
somewhat longer than 13 subunit s (some results sugges t 
a structure closer to a 28/13 helix ijEgelman et al.l [l982'0 
could also partially explain the deviation. 

An issue that has been much discussed is the contri- 
bution of Brownian motion and the power stroke to the 
total step size. With the geometric data used in this 
study, the power stroke, i.e., the distance of the lever 
arm tip movement between the states ADP.Pi and ADP, 
is about 31nm, or 5mn less than the average step size. 
Note that the second, smaller power stroke connected 
with ADP release does not contribute to the step size 
because it is normally followed by the detachment of the 
same head. Its function could be suppressing premature 
dissociation before the lead head binds and thus improv- 
ing the processivity. The remaining 5nm can be overcome 
by Brownian motion before the lead head binds. How- 
ever, at low loads, the binding of the lead head does not 
move the load, but rather stores the energy into bent 
lever arms. This energy gets released when the rear head 
detaches, which leads to an elastic power stroke imme- 
diately preceding the power stroke upon Pi release. At 
higher loads the situation is different, because the 5nm 
load movement occurs when the lead head binds. In nei- 
ther case we expect the 5nm power stroke to be resolvable 
under normal conditions because it always immediately 
precedes or follows the large power stroke. However, it 
is possible that the substeps become observable in the 
presence of chemical s that slow down the power stroke 
i|Uemura et"alll2004 . 

In order to fully re produce the substeps as reported by 
lUemura et alJ l)2004l) . some modifications would be nec- 
essary to the model. First, part of the power-stroke would 
have to occur immediately upon Pi release, resulting in a 
lever arm move of about 12 nm (first substep). This step 
would need a very strong force-dependence in its transi- 
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tion rate (activation point near the final state). The sub- 
sequent longer power stroke (ADP'^ADP) would then 
need a slower rate (~ 200 s -1 ) with less force dependence 
(activation point close to the initial state) . However, the 
finding that the substep position is independent of force 
remains difficult to explain, because the substep involves 
transition between a stiff configuration, bound on both 
heads, and a more compliant state, bound on a single 
head. 

The main value of both models (4- and 5-state) is that 
they provide a quantitative explanation of the coordi- 
nated head-over-head motility of the dimeric molecule 
while using only the properties of a single head as in- 
put. Both models also explain the observed force- velocity 
curves at high and low ATP concentration and the ef- 
fect of additional ADP, but these features already reveal 
some testable differences between the two scenarios. One 
of them is the shape of the force-velocity curve. In the 

4- state scenario the reverse power-stroke needs the re- 
binding of a phosphate molecule. This makes the cutoff 
behavior at high loads dependent on the Pi concentra- 
tion: the velocity drop is more gradual at low, but might 
become sharper at high Pi concentrations (Fig. |Hf3). In 
the 5-state scenario the velocity decline is more abrupt 
regardless of the Pi concentration. This is the first sug- 
gestion how experiments with improved precision and a 
wider range of chemical conditions could help distinguish- 
ing between the two scenarios. 

The main difference between the two scenarios is the 
predicted shape of the run length. Because the disso- 
ciation can take place on three different pathways, its 
rate depends on a number of parameters, of which a few 
cannot yet be determined by other methods. In the 4- 
state model the dissociation rate at high loads is domi- 
nated by detachment of a head in the ADP.Pi state and 
it therefore depends on the ratio k°_ A /k°_ Pi (Eq. I14fl . A 
strong increase with the load is characteristic for the 4- 
state model, because the load slows down the phosphate 
release and prolongs the dwell time in the state that is 
most vulnerable to dissociation. Dissociation at negative 
(forward) loads is dominated by pathways 2 (ATP medi- 
ated actin release in one head before the other head has 
bound) and 3 (dissociation of a head with ADP). In the 

5- state model all three pathways can contribute towards 
the dissociation rate, but there is no significant increase 
for positive loads - in fact, the dissociation rate can even 
decrease. 

The run length shortens with an increasing ADP con- 
centration in both scenarios. The decrease in run length 
is weaker than the decrease in the veloci ty (Fig . I12L 
which is consistent with recent observations ijBaker et all 
120041) . However, we cannot reproduce the reported com- 
plete saturation of run le ngth at hight ADP concentra- 
tions. iBaker et al. (2004) explain this saturation with a 
big difference (50-fold) between the attachment rates of 
the lead head depending whether the trail head is in the 
ADP or apo state, which we currently cannot reproduce 
with the relatively small power stroke (10°) upon ADP 



release in our model. 

An interesting difference between the 4- and the 5-state 
model is also that the 5-state model allows backward 
steps at high loads (above the stall force), while the 4- 
state model predicts rapid dissociation. In general, there 
are three possibilities how backward steps can occur: (i) 
The motor hydrolyzes ATP, but runs backwards, (ii) 
The motor slips backwards without hydrolyzing ATP — 
this is the case in our model, (iii) The motor synthesizes 
ATP from ADP and phosphate while being pulled back- 
wards, as ass umed by tightly coupled stoch astic stepper 
models (e.g.. iKolomeiskv and Fishier! |2003|) . It is possi- 
ble to test these three possibilities experimentally: If (i) 
is the case, the backward sliding velocity should show a 
Michealis-Menten type dependence on ATP concentra- 
tion. This mechanism would, however, require an even 
looser mechanochemical coupling, so that not only the 
release of Pi, but also the release of ADP and binding 
of ATP would be possible without completing the power 
stroke. In case (iii) it should depend on ADP as well as 
on Pi concentration, but not on ATP. In case (ii), which 
is favored by our study, the backward stepping occurs 
when both heads have ADP bound on them and they 
successively release actin at the lead position and rebind 
it at the new trail position. Even though this stepping 
requires no net reaction between the nucleotides, a cer- 
tain (low) ADP concentration is still required to prevent 
the heads from staying locked in the rigor (no nucleotide) 
state. 

The application of the elastic lever-arm approach de- 
veloped here should not be limited to simple geome- 
tries and longitudinal loads. A natural extension of the 
present work will be the influence of perpendicular forces 
on the activity of the motor. One will also be able to 
study the stepping behavior in more complex geometries, 
for example when passing a branching si t e indu ced by the 
Arp2/3 complex ijMacheskv and Gouldlll999|) . 

After completion of this manuscript, i t has been 
brought to my attention that lLan and Sul Eiol have 
also published a model for myosin V, based on the elas- 
ticity of the lever arm. In contrast to our model, they 
do not describe it as an isotropic rod, but use a weaker 
in-plane stiffness, combined with a strong (phenomeno- 
logical) azimuthal term that prevents binding of both 
heads to adjacent sites on actin. Another difference is 
that their study explicitly excludes dissociation events, 
whereas we use the dissociation rate to determine some 
of the model parameters. 
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Appendix 

Numerical solution for the lever arm shape 

The aim of this calculation is to determine the shape 
of the dimeric molecule for a given set of binding sites 
(trailing head bound on the site with the index i±, lead- 
ing head with 12), nucleotide states, which determine the 
lever arm starting angles cj>i and 4>2 1 and a given external 
load F. 

We start this task by deriving a function that numer- 
ically determines the endpoint of a lever arm as a func- 
tion of the force acting on it: Xj(Fj,(f>j) (j = 1,2). The 
shape of the whole molecule can then be determined nu- 
merically from the conditions that the endpoints of the 
two lever arms coincide, Xi = X2, and from the force 
equilibrium in that point 



F 2 = -Fe x 



(15) 



In many cases the function Xj will have more than one 
solution. Then we solve the system with all possible com- 
binations and then choose the solution with the lowest 
energy U = U\ + U2 + Fx, where U\ and U2 denote the 
energy stored in the distortion of each lever arm and Fx 
the work performed against the applied load. 

For a head bound at site i, the position of the proximal 
end of its lever arm in Cartesian coordinates reads 



and its initial tangent 



ia + S 
-iMn(0) 
i?cos(6>) 



cos(0) 
— sin(</>) sin(0) 
sin(</>) cos(0) 



(16) 



(17) 



where <p is the lever arm tilt (a function of the nucleotide 
state), S is the relative position of the lever arm proximal 
end (0 or 3.5 nm) and 8 is the azimuthal angle of the 
actin subunit to which the head is bound, 9 = 9oi with 
9q s» yjj x 360° s» 166°. The helix rise per subunit is 
a = 2.75 nm. 

If the force F acts on a lever arm that leaves the head 
in the direction t°, the whole lever arm will be bent in a 
plane spanned by the vectors t° and F. We can introduce 
a new two-dimensional orthogonal coordinate system in 
this plane, so that 



F,. 




F = I V j (18) 

Fy=Ft F x = \F-io{Fi )\ (19) 

In this coordinate system the shape can be determined 




FIG. 15 The shapes of an elastic beam anchored at one end 
and pulled by a given force F on its other end. The dashed 
line shows the unloaded beam. According to the sign of the 
initial curvature and the final angle (j>L the solutions can be 
divided into 4 classes. The beam corresponds to the myosin 
V lever arm, which is anchored in the head at one end and 
connected to a flexible joint at the other end. Note that the 
bending shown is exaggerated in comparison with realistic 
dimer configurations. 



by solving the equations 



M(s) = F A (x(L) - x(s)) = EI 



d(f>(s) 
ds 



dx ~ 
ds 



l = ( sin(» 
1 cos(0) 



(20) 
(21) 



with the boundary condition 0(0) = 0. The symbol "A" 
denotes the outer product, which is the out-of-plane com- 
ponent of the vector product. If we differentiate Eg. 1201 
by 4> we get 



El- 



ds 2 



-F x cos(0) + F y sin(0) 



(22) 



Through partial integration and taking into account the 
boundary condition M(L) = 0, we finally obtain 

EI 

IT \d S/ 

— F x (sm(/)L — sin0) + F y (cos (f>L — cos</>) (23) 

>L - " \ . I , 



= F sin 



sin 



Here we introduced the force angle <j>p, so that F x = 
F sm((f>F ) and F y — F cos(4>p). 

Because of the ambiguity of a quadratic equation, 
Eq. [23 generally has two solutions for a given set of 
values for <p(s), F, 4>l and <pp- As we have defined 
the coordinate system in a way that F x > 0, we have 
< <I>f < 7T- We also restrict ourselves to solutions 
with |0(s)| < 27r, i.e., we do not consider any spiraling 
solutions, because they always have a higher bending en- 
ergy than the straighter solution with the same endpoint. 
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There are four classes of functions cj>(s) that satisfy the 
condition that the RHS of Eq. [231 be positive: 



Solution <j>L 4>(s — ► 0) 



conditions 



I 

II 
III 

IV 



< 



< 



b F - 2tt < (j) L < 2(j) F - 2tt 

< < 

b. F - 2tt < 4> L < 24> F - 2tt 



The solutions III and IV have a turning point at </>o = 
— 2(7r — 4>f) — 4>l, where d(j)/ds changes sign. Ea. l2~3l can 
finally be transformed to 



L 
L 




I(<p L ) (cases I and II) (24) 

(2/(0 o ) + J(0 L )) (cases III and IV) 

(25) 

-1/2 



Note that for classes II and III the RHS of Eq. EH is 
not monotonous in </>l and there can be two solutions for 
a given L. Taking this into account, we obtain a total 
of up to 6 solutions. A situation in which all cases are 
represented is shown in Fig. 1151 

The configuration of the dimer is determined by solving 
Eq. El for all possible combinations of modes and taking 
the one with the lowest potential. The numerical inte- 
gration and solution were performed using NAG libraries 
(Numerical Algorithms Group) and the 3-d graphical rep- 
resentation of the calculated shapes was made with POV- 
Ray (www.povray.org). 
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